function dt_stumps = make_dt_stumps(X, Y, depth_limit, n_features)

feature_comb = combnk(1:size(X,2), n_features); 

temp = randperm(size(feature_comb, 1)); 
feature_comb = feature_comb(temp(1:30), :); 
dt_stumps = cell(size(feature_comb,1), 1); 

for i = 1:size(feature_comb, 1)
    fprintf('%d:', i); 
    dt_stumps{i} = boost_dt_train_multi(X(:,feature_comb(i,:)), Y, depth_limit);
end 

% dt_stumps = zeros(size(X,1), 10, n_features); 
% 
% for i = 1:n_features
%     garbage = rand(size(X,1), 10); 
%     dt_stumps(:,:,i) = bsxfun(@rdivide, garbage, sum(garbage, 2)); 
% end 